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SUMMARY 


INTRODUCTION 

In  modern  airborne  data  prooesalng.  It  la  imperative  tbat  the  shortest 
possible  solution  for  mathematical  problems  be  presented  to  the  computer. 

It  has  been  observed  that  the  mathematics  Involved  In  a  "stral^tforward” 
approach  to  some  AShT  geometrical  problems  would  strain  the  capabilities 
pf  a  small  scale  digital  computer.  To  alleviate  this  situation,  a  coordi¬ 
nate  transformation  methodology  has  been  developed  which  shows  great 
promise.  !Ihe  method  for  solving  problems  pertinent  to  ASW  presented  In 
this  report  Is  based  on  a  coordinate  transformation  technique.  A  proce¬ 
dure  employing  these  mathematical  considerations  is  developed  in  appen¬ 
dices  A,  6,  and  C. 

Two  applications  are  presented  as  a  testimonial  for  the  utility  of  the 
transformation  technique.  Application  I  is  that  of  determining  the  coordi¬ 
nates  of  Intersection  of  an  ellipse  and  a  confocal  hyperbola.  Application 
II  is  a  similar  determination  for  the  Intersection  of  two  ellipses  having 
one  common  focus. 

SUMMARY  OF  RESULTS 

1.  The  prescribed  method  effects  an  efficient  solution  (in  so  far  as 
digital  computation  is  concerned)  for  the  situations  exemplified  by  appli¬ 
cations  I  and  II.  !Ihe  validity  of  the  equations,  associated  with  these 
applications,  has  been  ascertained  by  programming  them  for  various  com¬ 
binations  of  their  pertinent  parameters. 

2.  It  Is  implicit  in  appendix  A  and  explicit  in  appendices  B  and  C  that 
even  though  trigonometric  and  hyperbolic  functionr  are  inherent  to  the 
development  of  the  proposed  technique  they  need  never  be  contended  with 
directly. 

3.  VSien  digital  computation  Is  to  be  considered,  the  transformation 
technique  Is  preferable  to  the  cosine  law  technique  for  the  solution  of 
the  problem  posed  by  application  I  ( see  appendix  D) . 

U.  Ihe  transformation  technique.  Illustrated  in  application  II  Is  pref¬ 
erable  to  a  competitive  method  (ellipse-ellipse  solution,  see  appendix  D) 
idien  digital  computation  Is  concerned. 

CONCLUSIONS 

1.  The  proposed  transformation  technique  leads  to  a  more  efficient  method 
of  solution  (with  regard  to  digital  computation)  for  the  problem  posed  in 
application  II. 
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2.  Although  a  cursory  comparison  of  the  transformation  technique  and 
the  cosine  law  technique  (application  I)  may  suggest  a  slight  preference 
for  the  latter,  on  the  grounds  of  general  familiarity,  the  programming 
analysis  of  appendix  D  showed  the  former  to  be  superior  vhere  digital 
eonputation  is  concerned.  (It  should  be  noted  that  a  digital  computer 
couldn't  care  less  about  the  intrinsic  nature  of  the  philosophy  pertinent 
to  a  method  of  solution.) 

3.  Appendix  D  shows  that  the  transformation  technique  (application  II) 
is  preferable  to  a  competitive  method  where  digital  computation  is 
concerned. 

ii.  Applications  I  and  II  demonstrate  the  simplifications  that  may  be 
expected  from  the  utilization  of  transformation  techniques  for  the  solution 
of  simultaneous  second-degree  equations  representing  intersecting  conic 
sections  which  arise  in  the  data  processing  of  many  tactical  situations. 

It  is  expected  that  a  general  solution  employing  such  techniques  might 
lead  to  the  optimum  means  of  programming  an  airborne  digital  computer  to 
solve  problems  of  this  nature. 

RECQMMBNBATIONS 

1.  It  is  recommended  that  the  method  illustrated  by  application  II  be 
employed  for  the  digital  solution  of  the  pertinent  tactical  situation. 

2.  It  is  recommended  that  the  transformation  technique  be  considered 
for  application  to  other  tactical  situations. 

3.  It  is  recommended  that  other  transformation  techniques  be  studied 
for  the  purpose  of  reducing,  still  further,  the  amount  of  work  to  be 
done  by  a  digital  computer. 


-  ii  - 


NADC-AW-6239 


TABLE  OF  CONTENTS 

Page 

SUMMARY .  i 

Introduction .  i 

Summary  of  Results .  1 

Conclusions .  1 

Recommendations .  11 

APPLICATIONS  .  1 

Application  1 .  1 

Application  II  .  2 

Solution  Summary . «...  3 

Figures 

1  Geometry  Pertinent  to  Application  I  .  1 

2  Problem  Geometry  (Application  II)  .  .  .  . .  2 

APPENDICES 

A  Ihe  Coordinate  Transformation  Technique .  A-1 

B  Mathematics  Associated  with  Distance  Measuring 

Equipment  (DHE) . B-1 

C  Mathematics  of  Transformation  Technique  Applied  to 

Intersecting  Ellipses .  C-1 

D  Programming  Analysis  of  Applications  I  and  II  .  .  .  D-1 


-  Ill  - 


IIADC-AW-6239 


APPLICATIONS 


Applloctlon  I  Illustrates  a  method  (based  on  a  transformation  tech¬ 
nique)  for  the  determination  of  a  target  position  relative  to  a  coordi¬ 
nate  frame  defined  by  two  fixed  points  when  the  distance  between  these 
points  and  the  undirected  distances  from  each  point  to  the  target  are 
known.  Die  final  equations  shown  herein  are  derived  In  appendix  B. 

Application  II  Illustrates  a  method  for  the  determination  of  the 
points  of  Intersection  of  two  ellipses  with  a  common  focus.  (Ihese 
Intersections  may  also  be  defined  by  the  Intersection  of  one  of  the  above 
mentioned  ellipses  and  one  branch  of  a  hyperbola. )  Ihe  equations  employed 
herein  are  derived  In  appendix  C. 

APPUCATION  I 


Problem I  To  determine  the  coordinates  (u,  v)  of  a  point  P  when  the  ranges 
(RX  and  R2)  to  two  fixed  points  are  specified.  Figure  1  Illustrates  the  ■ 
pertinent  geometry. 


% 


Pi 

P2 


u 


fixed  point  1 
fixed  point  2 


Line  through  Px  and  P2  defines  u  axis. 

Perpendicular  bisector  of  line  Joining  Px  and  P2  defines  the  v  axis. 


FIOURE  1  -  Qeometry  Pertinent  to  Application  I 


Solution  Suimiary 
u  -  (Ajj)  (A,)  Y 

^  -  ^(Bh)  (Be)  ^ 

where  Ag  -  Bg  •  ±  .^A*g  -  1 

R3 

Ah  -  Bh  -  Vl 
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axul  the  x^eolutlon  of  the  enaulng  ambiguity  Is  oonaldered  to  be  a  separate 
problem}  however,  In  an  applloatlon  auoh  aa  IXfE  a  prior  knowledge  of  the 
approxlMte  location  of  the  point  P  will  auffloa  to  yield  the  deaired 
result. 

APPLICJATION  II 

Problemt  To  determine  the  location  of  a  target  relative  to  two  signal 
defection  polnta  P]^  and  P2. 

QLveni  1,  The  coordinates  of  two  signal  detection  points  and  that  of 
the  signal  source. 

2.  TWO  numerical  values:  and  k2. 

where  kj,  ••  difference  between  the  time  of  arrival  of  the  direct  and 
(1  ■  1  or  2)  reflected  signals  (at  signal  detection  point  l)  times 
the  velocity  of  the  signal  in  the  applicable  medium. 


Assumption:  Rectilinear  signal  propogatlon. 

Figure  2  Illustrates  the  pertinent  problem  geometry. 


V 


FIGURE  2  -  Problem  Geometry  (Application  II) 

Given:  a  k^ 
b.  k2 
c 

Find:  uj 
▼T 
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Solution  0— Mcy 
OoiVUtM 

•  -  ^(b  *•/  *  <f 

» -  y(b  -.)•*<<• 

oi-V^ 

4.!!^ 

^  a 
b' 

a 

o‘  -  2 

a 

I  I 

,  -  «J2 

•h-— r“ 

R  -  2  [b'  (d[  .  dj)  -  d’  -  d^j 

P  ■  2c^  j^d^  “  ^*2] 

M  -  (b-)*  +  (c')*  -  1 

Q  -  d^  [m  ♦  (dj;)*]  -  di  [m  +  (d^)"] 

(b^)»  -  1  -  (ajj)- 
A  -  (bjj)*  -  H?  (ajj)* 


B  ■  -2  QRajJ 
C  -  ~P>  (bjJ)*  -  (? 
then,  for  A^O 
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'"1,2  "  'H 


-B  ±  t  -  U  AO 


2A 


'11,2  -  •  'll, 2 


'^1,2  "  ® 


I 


If  A  -  0 


then 


ui,2  ■  0 

-H  ±  -/rf*  -  li  (E- 

n,2  "  - ^ 

2Q 


where 


0  -  li  (c  -  dj^) 


•  adj^  or  ♦  a 

M*l?  +  0?*  -  <i*-dj_ 

H  -  -U  e  M 

Z  •  rf  -  lidfa* 

Reenltt 

Ambiguous  solution  It 
Ambiguous  solution  2i  U2,  ^2 

VIhere  the  resolution  of  the  ensuing  ambiguity  Is  considered  to  be 
separate  problem. 

*  The  subscript  1  Implies  the  use  of  the  {*)  sign. 

Ihe  subscript  2  Implies  the  use  of  the  (-)  sign. 
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THE  COORDINATE  TRANSFORMATION 

TECHNIQUE 

This  appendix  presents  the  use  of  a  transformation  teclmique  to 
effect  a  more  efficient  solution  for  a  particular  problem.  The  first 
part  discusses  transformations  in  general  and  develops  those  particular 
transformation  equations  which  are  employed  thTOu^out  this  reportj  the 
latter  part  discusses  the  transformation  technique  ^diich  is  employed  for 
the  determination  of  the  points  of  intersection  of  two  hyperbolas.  (Since 
this  problem  is  of  a  complex  nature  and  presents  no  foreseeable  applica¬ 
tion,  the  reader  may  omit  it  without  lose  of  continuity. ) 

The  solution  of  many  problems  becomes  very  much  simpler  idien  the  for¬ 
mulas  are  expressed  in  terms  of  a  suitable  coordinate  system.  Many  tac¬ 
tical  situations  involve  the  solution  of  simultaneous  second-degree  conic 
equations  such  as  hyperbolas.  The  object  of  this  report  is  to  investi¬ 
gate  the  feasibility  of  obtaining  a  more  convenient  coordinate  transfor¬ 
mation  for  these  systems,  l.e.,  reduce  the  problem  to  the  simultaneous 
solution  of  linear  equations  and  make  lighter  the  task  of  an  airborne 
digital  computer. 

MATHEMATICAL  BACKfflOUND  (COORDINATE  TRANSFORMATIONS) 

Consider  two  planes,  one  in  idilch  a  point  is  specified  by  the  numbers 
X,  y  and  the  other  in  lAich  a  point  is  specified  ly  the  numbers  u,  v.  A 
transformation  sets  up  a  correspondence  between  points  in  the  x,  y  plane 
lying  in  a  certain  domain  and  points  in  a  certain  domain  in  the  u,v  plane. 
This  correspondence  may  be  specified  by  the  transformation  equations: 

u  “  f  (x,  y),  V  =  g  (x,  y) 

^ere  the  functions  f  and  g  are  required  to  be  single  valued  so  that  the 
point  (u,  v)  is  uniquely  determined  when  a  point  (x,  y)  is  given.  It  is 
convenient  to  call  one  of  the  planes,  say  the  x,  y  plane,  the  original 
plane,  and  the  u,  v  plane  the  image  plane,  and  say  that  the  point  (u,  v) 
is  the  image  of  the  point  (x,  y).  If  the  functions  f  and  g  are  single 
valued,  there  is  only  one  image  point  for  every  point  in  the  original 
plane.  This  does  not  preclude  the  possibility  that  the  point  (u,  v)  may 
not  be  the  image  of  several  points  in  the  x,  y  plane.  If,  however,  the 
point  (u,  v)  is  the  image  of  only  one  point  in  the  original  plane,  the 
transformation  is  said  to  be  reversible.  The  condition  for  this  is  that 
the  functions  f  and  g  possess  unique  inverses.  The  condition  for  this  is 
simply  that  the  o^acobian  of  the  transformation  not  vanish  over  the  domain 
under  consideration.  Tnus: 
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a  u 

a  X 

a  V 
a  X 


a  u 
a  y 

a  V 
a  y 


TECHNIQUE 

A  procedure  is  developed,  based  on  conformal  mapping  procedures,  for 
rapidly  determining  the  point  of  intersection  of  a  pair  of  hyperbolas. 

![he  method  to  be  described  is  the  geometrical  analogue  of  operational 
calculus  techniques  such  as  the  Laplace  transformation  vdiere  dlfferen/blal 
equations  in  time  are  transformed  into  algebraic  equations  in  terms  of  a 
secondary  variable.  This  simpler  problem  is  solved  in  terms  of  the  sec¬ 
ondary  variable  and  then  transformed  back  into  the  original  independent 
variable  to  obtain  the  solution  of  the  original  differential  equation. 


Ihe  prescribed  procedure  calls  for  the  mapping  of  conic  sections  in 
the  u,  V  plane  into  stral^t  lines  in  the  x,  y  plane,  determining  the 
intersecting  point  of  these  strai^t  lines  in  the  x,  y  plane  and  then 
mapping  this  point  back  into  the  u,  v  plane  for  the  desired  solution. 


DEVELOPMENT  OF  EQUATIONS 

Uie  transformation  W  =  sin  Z  will  be  employed  for  this  particular 
study. 


W  «  u  +  Iv  -  sin  Z 

-  sin  (x  +  iy) 

“  sin  X  cos  iy  +  cos  x  sin  iy 

=  sin  X  cosh  y  +  cos  x  [  i  sinh  y] 

=  sin  X  cosh  y  +  1  cos  x  sinh  y 

Therefore,  the  transformation  W  »  sin  Z  can  be  written 


u  *  sin  X  cosh  y  (A-l) 

V  =  cos  X  sinh  y  (A-2) 

and  the  line  x  =  k  =  constant,  >diere  -^  <  k  <  ^,  maps  into  the  curve 
idiose  parametric  equations  are 

u  ■  sin  k  cosh  y  (A-3) 

V  “  cos  k  sinh  y  (A-]|) 
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As  an  example,  if  k  =  30*  =  0.52360  radians  then 
u  »  0.50000  cosh  y  =  f]_(y) 

V  *  0.86603  Binh  y  f2(y) 

Table  A-I  lists  values  of  u  and  v  for  the  various  values  of  y  in  these 
equations . 

TAB.  LE  A-I 
u  =  fl(y)»  V  =  f2(y) 

£  cosh  y  u  °  0.50000  cosh  y  sinh  y  v  ■  0.86603  sinh  y 


0.0 

1.00000 

0.50000 

0.00000 

0.00000 

0.2 

1.00201 

o.5iooli 

0.2013li 

0.17593 

O.ii 

1.08107 

0.5ii05ii 

0.1*1075 

0.35572 

0.5 

1.12763 

0.56382 

0.52110 

0.45130 

0.7 

1.25517 

0.62759 

0.75858 

0.65695 

0.9 

1. 1^3309 

0.71655 

1.02652 

0.88900 

1.1 

1.66852 

O.83U26 

1.33565 

1.15672 

1.3 

1.97091 

O.985I16 

1.69838 

1.1*7085 

Using  the  identities 
sinh  (-y)  =  -  sinh  y 
cosh  (-y)  =  cosh  y 
It  can  be  seen  that 
u  =  0.50000  cosh  y 

V  =  -0. 86603  sinh  y 

for  negative  value  of  yj  therefore,  the  curve  defined  by  the  parametric 
equations  (A-3)  and  (A-U)  will  be  symmetric  with  respect  to  the  u  axis. 

The  curve  for  x  =  +30®  is  shown  in  figure  A-1. 

Likewise,  if  k  =  -30"  =  -0.52360  radians  then 

u  =  sin  (-30)  cosh  y 

V  *  COB  (-30)  sinh  y 
and  using  the  identities 

-  A-k  - 
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sin  (-x)  =  -  sin  x 

cos  (-x)  =  cos  X 
then, 

u  *  -  sin  30*cosh  y 

V  =  cos  30*sinh  y. 

Table  A-I  can  be  employed  to  compute  this  curve  which  is  also  shown 
in  figure  A-1. 

To  demonstrate  mathematically  that  parametric  equations  (A-3)  and 
(A-li)  represent  a  hyperbola,  square  (A-3)  and  (A-lj.)  respectively  to 
obtain 


if*  =  sir?  k  cosh?  y 
1?  =  co^  k  sinl^  y. 
Substitute  the  identity 
cosl^  y  =  1  +  sinh?  y 
into  equation  (A-5)  to  obtain 
i?  =  sin®  k  (1  +  sinh®  y). 


(A-5) 

(A-6) 

(A-7) 

(A-8) 


Now,  substitute  the  expression  for  sinh*  y  from  equation  (A-6)  into 
equation  (A-8)  to  obtain 

V® 

t?  =  sir?  k  (1  +  — - — )  (A-9) 

cor  k 


idiich  reduces  to  the  standard  form  of  the  hyperbolic  equation  as  follows: 
v?  V®  i;?  V® 

— - - - —  “  1  “  T  ■  T  (A-10) 

sir?  k  co^  k  a®  t? 


letting  sin  k  -  a 
cos  k  -  b. 

For  this  transformation  it  has  been  demonstrated  that  the  following 
conditions  apply: 

If  k  >  0,  the  vertical  line  in  the  xy  plane  maps  into  the  ri^t-hand 
half  of  a  hyperbola  in  the  uv  plane; 

If  k  <  0,  the  vertical  line  in  the  xy  plane  maps  into  the  left-hand 
half  of  the  hyperbola  in  the  uv  plane; 
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-  Plot  of  Paraaetrio  Equations 
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The  mapping  is  one  to,  one.  The  points  w  =  ±1  4re  the  foci  of  the 
hyperbola. 

From  the  above  it  can  be  seen  that  vertical  lines  in  the  xy  plane 
map  into  hyperbolas  in  the  uv  plane  under  the  transformation  w  ■  sin  Z. 

Conversely,  hyperbolas  in  the  uv  plane  map  into  vertical  lines  in  the 
xy  plane,  since,  if  w  is  a  function  of  Z,  it  follows  inversely  that  Z  is 
a  function  of  w. 

E3CAMPLE  PROBLEM  -  INTERSECTING  HYPERBOUS 

An  example  of  the  use  of  the  w  >  sin  Z  transformation  in  solving  prob¬ 
lems  involving  intersecting  hyperbolas  follows; 

The  problem  which  will  be  solved  here  is  to  determine  the  point  of 
intersection  of  hyperbola  branches  No.  1  and  No.  2  as  illustrated  in 
figure  A-2. 

As  can  be  seen  from  figure  A-2,  hyperbola  branch  No.  1  may  be  most 
easily  defined  relative  to  the  uv  coordinate  frame  with  points  f  and  f ' 
as  foci,  and  similarly  hyperbola  branch  No.  2  may  most  easily  be  defined 
relative  to  the  u'v'  coordinate  system  with  points  f  and  f ' '  as  foci. 

Also,  for  this  problem  the  u'v'  coordinate  system  has  been  translated 
and  rotated  throu^  an  angle  of  90  degrees  relative  to  the  uv  coordinate 
system. 

The  problem  has  also  been  designed  so  that  the  distances 

|0f  I  *  |of'  I  =  |0'f  I  =  jO'f"  I  =  1. 

Problem  Particulars 


The  branch  of  hyperbola  No.  1  plotted  in  figure  A-3  from  data  pre¬ 
sented  in  table  A-II,  with  foci  (1,  0)  and  (-1,  0)  in  the  uv  plane,  has 
been  designed  to  pass  through  the  point  (3,  2) .  The  hyperbola  Is  defined 
as  the  locus  of  a  point  idilch  moves  so  that  the  difference  of  its  undi¬ 
rected  distances  from  two  fixed  points  is  a  constant,  the  value  of  idiich 
is  2a.  The  value  of  a  is  shown  in  figure  A-3. 

Choosing  the  axes  as  shown  in  figure  A-3  and  denoting  by  2a  the 
absolute  value  of  the  difference  between  the  distances  f'P  and  fP  it 
follows  that 

-  (-c))*  +  V®  -  ^{u  -  cf  +  V*  =  2a 
^3  -  (-1)  *•  +  2®  -  y(3  -  !)*  +  ^  =  2a 
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FIOORE  A-3  -  Set-Up  for  Ifyperbola  No.  1 
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-^16  +  li  -  1;  +  li  =  2a 

^  -  li.U72  -  2.828  -  -  2a 

a  =  0.822. 

Also,  from  the  relationship 
c®  =  +  a® 

it  follows  that 

t?  =  1  -  (0.822)®  =  1  -  0.676 

^  =  0.32li, 
and  from 

^  -  1 

a®  ’  I?  ’ 

the  equation 


=  b»  (-  -  1)  =  0.321;  -  1) 

a®  0.676 

is  derived  and  employed  to  develop  table  A-II  from  >rtiich  figure  A-3  has 
been  plotted. 

TABLE  A-II 

Points  which  define  the  locus  of  hyperbola 
No.  1  in  the  uv  plane  shown  in  figure  A-2 


u 

V 

0.822 

0.000 

1.000 

±0.395 

2.000 

±1.263 

3.000 

±1.999 

Similarly,  as  can  be  seen  from  figure  A-2,  hyperbola  No.  2  is  defined 
with  the  points  f  and  f>'  as  fooi.  Ihls  hyperbola  is  set  up  so  that  it 
intersects  hyperbola  No.  1  at  one  of  the  points  shown  in  table  A-II  such  as 
(2,  1.263)  in  the  uv  plane.  The  most  convenient  means  of  obtaining  an 
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equation  for  this  hyperbola  vould  be  to  establish  the  u'v'  coordinate 
system  as  shown  In  figure  A-2,  and  derive  the  mathematical  expression 
for  hyperbola  No.  2  In  terms  of  the  u'v'  coordinate  system. 

AlsOj  the  point  (2,  1.263)  In  the  uv  coordinate  system  Is  the  point 
(0.263j  -l)  In  the  u'v'  coordinate  system.  The  set-up  for  hyperbola  No. 

2  relative  to  the  u'v'  system  Is  shown  In  figure  A-k  as  plotted  from  data 
developed  for  table  A-III  as  follows: 

Again,  from  the  definition  of  a  hyperbola 

^(u'  +  c)®  +  v'*  -  ^ (VL*  -  c')®  +  v'  »  ±2^^' 

^[0.263  -  (-1)]*  +  1®  -  .yd  -  0.263)*  +  1®  -  2a' 


then 


a'  =  0.181* 

and  from  the  relationship 
cl*  -131*+  a'® 

It  follows  that 
ti*  -  c'®  -  a'® 

-  1  -  (0.181*)® 

=  0.966. 


Then,  from 


u' 


V'- 


a'*  b'® 

the  equation 


V'®  -  b'®  -  1)  -  0.966 

a'® 


u'- 


(0.181*)® 


Is  employed  to  develop  table  A-III. 
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TABLE  A-III 

Points  vdilch  define  locus  of  hyperbola  No.  2  in  u'v'  plane 


u*  V* 

0.l8ii  0.000 
0.263  ±0.998 
1.000  ±5.236 

To  Bunmarize  the  problem  initially  Illustrated  in  figure  A-2,  figure 
A-5  has  been  prepared  showing  the  point  of  intersection  of  hyperbolas  No. 
1  and  No.  2. 

Conventional  Method  of  Solution 


The  "strai^tforward"  method  of  solving  for  the  point  of  intersection 
of  hyperbolas  No.  1  and  No.  2  (as  shown  in  figure  A-5)  in  the  uv  plane 
would  be  to  solve  simultaneously: 


—  »  1  ■  ^  V* 

I?  0.676  0.32li 


(A-11) 


111  *  1  = 

b'®  0.03ii0 


0.966 


(A-12) 


Biis  technique  requires  that  the  equation  for  hyperbola  No.  2  be 
expressed  in  terms  of  the  u  and  v  coordinate  axes.  This  is  effected,  in 
general,  by  both  a  rotation  and  a  translation  transformation  hereby  the 
general  form  for  the  above  equations  become  respectively. 


Ai  V?  +  Cl  +  Pi  =  0  (A-13) 

A2  +  B2  uv  +  C2  +  D2  u  +  E2  V  +  F2  =  0.  (A-lli) 


Ihe  prescribed  technique  for  the  solution  of  the  above  equations 
involves  considerable  effort  not  only  in  algebraically  solving  for  the 
desired  solution,  but  also  in  finding  the  proper  expression  for  the 
equations  relative  to  the  same  coordinate  i^stem. 

Ihe  remainder  of  this  appendix  Illustrates  the  proposed  conformal 
transformation  technique  for  the  rapid  determination  of  the  point  of 
interseotlon  as  shown  in  figure  A-5. 

DISCUSSION  OP  THE  'mANSPORMATlON  SOLUTION  TECHNIQUE 

Hyperbolas  No.  1  and  No.  2  of  figure  A-5  will  be  mapped  into  straight 
lines  in  the  ±y  plane.  Ihe  point  of  intersection  of  these  straight  lines 
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in  the  xy  plane  will  then  correspond  to  the  intersection  of  the  hyper¬ 
bolas  in  the  uv  plane^  and  the  inverse  transformation  of  this  point  back 
into  the  uv  plane  will  then  yield  the  desired  solution. 

Transformation  of  hyperbola  No.  1  from  the  uv  coordinate  system  to 
the  xy  coordinate  system  is  accomplished  as  follows! 

Equation  of  hyperbola  No.  1 

~  ■  1  »  ^  T* 

I?  sin®  k  co^  k 

by  letting  a  *  sin  k  and  b  =  cos  k  as  previously  shown  in  equation  (A-IO). 
Thus,  sin  k  =  a  =  0.822 

k  =  sin“^  0.822 

k  =  55*  17^ 

k  =  0.961i  radians 

Figure  A-6  is  a  pictorial  exposition  utilizing  the  above  computations. 

The  task  of  transf owning  hyperbola  No.  2  is  not  as  strai^tforward  as 
that  encountered  in  the  hyperbola  No.  1  case.  However,  vblle  not  gener¬ 
ally  applicable^  in  this  specific  case  the  problem  can  be  attacked  without 
effectively  rotating  or  translating  the  u'v*  coordinate  system  into  the 
uv  coordinate  system. 

The  equation  of  hyperbola  No.  2  relative  to  the  primed  system  may  be 
written  as  follows: 

u'®  ^  .  u»®  v»® 

a»®  b*®  six?  k'  co^  k' 

by  letting  a'  =  sin  k'  and  b'  =  cos  k'  as  previously  diown  in  equation 
(A-10) .  Thus, 

sir?  k'  =  (a®)  =  (0.l8i:)* 

sin  k'  =  O.I8I4. 

k'  =  sin"^  O.I8I: 

k'  =  10“  38' 

k'  -  0.186  radians. 
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Since  hyperbola  No.  1  curving  to  the  right  in  the  u,  v  plane  trans¬ 
formed  to  a  line  parallel  to  the  y  axis  in  the  x,  y  plane,  it  is  reason¬ 
able  to  expect  that  hyperbola  No.  2  curving  upward  in  the  u',  v'  plane 
will  transform  to  a  line  parallel  to  the  x  axis  in  the  x,  y  plane.  The 
transformed  equation  would  be  expected  to  be  as  follows  relative  to  the 
X,  y  plane 

y  -  constant  -  af  +  k‘. 

Where  o’  is  the  y  coordinate  value  of  the  intersection  of  the  x,  y 
equation  for  the  v'  axis  relative  to  the  u,  v  coordinate  frame  and  the 
line  No.  1  in  the  x,  y  plane.  The  transformation  equation  for  the  v' 
line  (v  »  l)  in  the  uv  coordinate  system  is  given  by  the  following 
equation: 

V  -  1  «  cos  X  slnh  y 


or 


-  =■  sec  X  “  sinh  y. 

cos  X 


Table  A-IV  lists  values  of  y  for  various  values  of  x  in  the  above 
equation. 


TABLE  A-IV 

X 


Radians 

Degrees 

sec  X 

1 

0.000 

0.0 

1.0000 

0.89 

0.262 

15 

1.0530 

0.91 

0.523 

30  ■ 

1.151:7 

0.99 

0.785 

1:5 

1.1:11:2 

l.ll: 

1.01:7 

60 

2.0000 

1.1:5 

1.309 

75 

3.8637 

2.06 

1.571 

90 

Table  A-IV  is  employed  to  plot  the  above  equation  in  figure  A-7  which 
is  an  illustration  of  how  both  hyperbolas  No.  1  and  No.  2  map  onto  the  xy 
plane. 

As  shown  in  figure  A-7,  a'  may  be  Interpreted  as  the  transformation 
of  a  from  the  uv  coordinate  system  into  the  xy  coordinate  system. 

The  point  of  Intersection  of  lines  No.  1  and  No.  2  in  the  xy  plane 
uniquely  corresponds  to  the  point  of  intersection  of  hyperbolas  No.  1  and 
No.  2  in  the  uv  plane. 
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The  xy  coordinates  of  the  point  of  intersection  in  the  xy  plane  are 
(0.96U3>  I.5156),  and  the  transformation  equations  introduced  as  equa¬ 
tions  A-1  and  A-2  may  be  employed  to  map  this  point  into  its  image  point 
in  the  UT  plane  which  in  this  case  is  the  desired  solution. 

u  =  sin  X  cosh  y  (A-l) 

=  sin  (O.96U3  rad)  cosh  I.516 

=  (0,822)  (2.395) 

=  1.969 

V  =  cos  X  sinh  y  (A-2) 

=  cos  (0.961i  rad)  sinh  (I.516) 

=  (0.569)  (2.177) 

*  1.2li0 

Therefore,  (1.969,  1.21^0)  is  the  computed  point  of  intersection  of 
the  two  hyperbolas  in  the  uv  plane. 

SUMMARY 

In  summary,  the  problem  of  finding  the  intersection  of  the  following 
two  hyperbolas  illustrated  in  figure  A-5  involves  the  following  steps. 
Given: 


Hyperbola  No.  1 

u? 

V* 

- =  1 

(Right  hand  branch  only) 

(0.822)* 

(0.570)* 

Hyperbola  No.  2 

u'* 

v* 

- =  1 

(U^per  branch  only) 

(0.181I)® 

(0.966)* 

Procedure: 

X  -  sin'^  0.822  =  55"  17' 
y  =  sin**^  (O.I8I;)  +  a' 

a'  =  sinh"^  - - -  =  1.330 

cos  (0.822) 

y  =  0.1856  +1.33  -  1.516 
u  =  sin  X  cosh  y  =  (0.822)  (2.395)  =  1.969 
V  =  cos  X  sinh  y  =  (0.570)  (2.177)  "  1.2l|0 
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APPENDIX  B 

MAffiEMATTCS  ASSOCIATED  WITH 
DISTANCE  MEASURING  EQUIPMENT  (l]ME)« 

(The  first  practicable  application  of  the  "Transformation  Technique"  is 
presented  herein^*-"-) 
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*  These  resxxlts  are  also  applicable  to  other  tactical  situations  having 
a  similar  geometrical  configuration. 

Authors  Note:  Since  the  simultaneous  solution  of  the  equations 
representing  two  intersecting  circles  yields  a  simpler  ablution  than 
those  discussed  herein,  this  appendix  is  presented  for  illustrative 
purposes  only. 
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MATHEMATICS  ASSOCIATED  WITH  DISTANCE 
MEASURING  EQUIPMENT  (DME) 

THE  NAVIGATIONAL  PROBLEM 

Precise  navigation  over  the  surface  of  the  ocean  is  difficult  to 
accomplish.  The  object  of  radio  DME  is  to  improve  the  navigational 
accuracy  of  aircraft  in  ASW  operations  by  electronically  locating  the 
aircraft  with  respect  to  sonobuoys  modified  to  operate  as  electronic 
marine  markers.  Radio  DME  systems  are  not  dependent  upon  good  weather 
conditions,  time  of  day,  the  operator's  eyes,  or  the  proximity  to  con¬ 
ventional  marine  markers. 

OBJECTIVE 

The  objective  of  this  appendix  is  to  determine  the  simpler  of  two 
mathematical  expressions  for  the  position  of  an  aircraft  relative  to  two 
sonobuoys . 

SYSTEM  GEOMETRY 

The  position  of  an  aircraft  may  be  uniquely  determined  relative  to  a 
coordinate  frame  defined  by  two  sonobuoys  as  follows j  a  strai^t  line 
throu^  the  sonobuoys  may  conveniently  be  taken  as  the  u-axisj  a  strai^t 
line  through  one  of  the  sonobuoys  and  perpendicular  to  this  u-axis  may  be 
taken  as  the  v-axis  as  shown  in  figure  B-1.  ^ 

PROBLEM  STATEMENT 

Two  simple  mathematical  expressions  for  determining  the  position  of 
an  aircraft,  relative  to  two  sonobuoys  may  be  derived  from  either  the 
"coordinate  transformation"  or  "cosine  law"  techniques. 

Coordinate  Transformation  Technique 

The  coordinate  transformation  technique  as  described  in  appendix 
A  may  be  employed  to  determine  the  aircraft  position  (P)  relative  to  the 
u-v  coordinate  frame.  This  technique  is  illustrated  in  figure  B-2, 
where  the  transformation  w  =  sin  z  (w  =  u  +  iv  and  z  =  x  +  iy)  maps  the 
upper  half  of  the  ellipse  from  the  u-v  plane  into  the  straight  line  E^p 
in  the  x-y  plane,  and  the  lower  half  of  the  ellipse  from  the  u-v  plane^ 
into  the  straight  line  Erg  in  the  xy  plane.  Also,  the  right  and  left 
halves  of  the  hyperbola  in  the  u-v  plane  are  mapped  by  this  transformation 
into  the  straight  lines  Hj^  and  H^,  respectively,  in  the  x-y  plane. 

Considering  the  symmetry  of  the  problem  and  taking  the  aircraft  to 
be  at  position  P  in  the  u-v  plane,  figure  B-3  may  be  constructed  to 
illustrate  the  application  of  the  transformation  technique  for  the 
determination  of  aircraft  position. 
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V 


Ihe  symbols  shown  abore  are  defined  as  follows! 

SB(l)  ■  Position  of  sonobuoy  1 
SB(2)  ■  Position  of  sonobuoy  2 

R;i^  >•  Distance  between  SB(l)  and  aircraft 

1 

R2  *  Distance  between  SB(2)  and  aircraft 
R^  >  Distance  between  sonobuoys 

P(uj  v)  -  Position  of  aircraft  relative  to  u,  y  coordinate  frame 

FIGURE  B-1  -  Aircraft  Position  Relative  to  a  Coordinate  Frame 
Established  by  Two  Sonobuoys 
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Point  P  is  unlquelj  napped  into  point  P' 


TLGOBE  B-3  -  Equations  tor  Detaminlng  Aircraft  Position  "P"  on  the  u-t  Coordinate  Fraae 


NADC-AW-6239 


In  appendix  A,  it  was  shown  that  the  coordinates  of  the  aircraft  in 
the  u-v  coordinate  frame  could  be  given  as  a  function  of  a  point  (x,  y) 

in  the  7tr7  coordinate  frame  where 
% 

u  »  sin  X  cosh  y 


and 


V  =  cos  X  sinh  y. 

Now,  it  is  possible  to  further  simplify  these  expressions  so  as  to 
avoid  the  need  for  evaluating  the  above  trigonometric  and  hj^erbolic 
functions . 

A  consideration  of  the  geometry  involved,  along  \d.th  the  equations 
given  in  appendix  A  will  reveal  that 

sin  X  =  Afj 

cosh  y  “  Ag 

cos  X  =  Bh 

sinh  y  *■  Bg 

where  the  subscripts  E  and  H  refer  to  the  defining  ellipse  and  hyperbola, 
respectively}  that  is^  2  A^  is  the  constant  referred  to  in  the  classical 
definition  of  the  hyperbola  shown  in  the  u,  v  plane  in  figure  B-2  vdiich 
passes  throu^  the  point  P. 

Noting  that  the  hyperbola  is  defined  as  the  locus  of  a  point  idiich 
moves  so  that  the  difference  of  its  undirected  distance  from  two  fixed 
points  is  a  constant,  the  problem  may  be  arranged  as  shown  In  figure 
B-3  with  the  scale  adjusted  so  that  R2/2  »  1,  as  required  by  the  theory 
given  In  appendix  A.  In  this  manner,  may  be  evaluated  as 

1/2(R|  -  Rp  and  Bjj  may  be  evaluated  as  (l  -  where 

,  Rl  I  ^2 

"1  ■  572  ■  572- 

Also,  2  Ag  is  the  constant  referred  to  in  the  definition  of  the  ellipse 
which  also  passes  through  the  point  P. 

Noting  that  the  ellipse  is  defined  as  the  locus  of  a  point  >dilch 
moves  so  that  the  sum  of  its  undirected  distances  from  two  fixed  points 
is  a  constant.  Again,  when  the  problem  is  arranged  as  shown  by  figure  ^3 
with  the  scale  adjusted  so  that  R^2  =  1,  Ag  may  be  evaluated  as 

1/2(R^  +  r') 
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and  Bp  may  be  evaluated  as 

where  and  Rg  are  given  above. 

Finally,  the  following  equations  may  be  utilized  to  find  the  aii*oraft 
position  "P"  in  the  u-v  coordinate  frame. 

u  -  (Ajj)  (Ag)  (R3/2) 

V  =  (Bh)  (Bg)  (R3/2) . 

Ihe  following  sample  problem  (figure  B-lt)  may  be  utilized  to  demon¬ 
strate  the  coordinate  transformation  technique. 

Sample  Problem 

Given: 

R^  -  10,000.000  yd 
Rg  =  8,08U.736  yd 
R3  -  2,000.000  yd 
Solution: 

Ag  •  |(10  +  Q.OQkiyS)  =  I  (18.08U736) 

=  9.042368 

Ajj  -  I  (10  -  8.084736)  =  I  (1.915264) 

=  0.957632 

^  (R3/2) 

-  (9.042368)  (0.957632)  (1,000) 

=  8,659.261 

®E  “  “v/  Ag  -  1 

*  -v/(9. 04236)*  -  1  =  -v/80. 78670096 
=  8.988142 
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If  -  10,000  yd 

R3/2  -  1,000  yd 

e  •  15*  00' 

Then,  u  -  cos  9  -  R3/2 

-  10,000  (0.96533)  -  1,000 

-  8,659.3 

and  ▼  ■  Rj^  Bln  15* 

-  10,000  (0.25882)  -  2,580.2. 

Also,  using  the  distance  formula 

“  V  (u  -  nyg)^  ♦  ^  -  y (7,659. 3)*  +  2,588.2)* 
“  65.3630  X  icf  -  8,081|.7. 

FIGURE  B-k  -  Sain)le  Problem  Arrangement  for  the 
Transformation  Technique 
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Br  =  y  1  -  A§ 

=  ^/l  -  (0.957632)*  1  -  0.9170590li7  =  .y^0829l095 

-  0.28799ii 

V  =  (Br)  (Br)  (R3/2) 

-  (6.98811^2)  (0.28799li)  (1,000) 

=  2,588.330 

Cosine  Law  Technique 


The  cosine  law  may  be  stated  as  follows:  "The  square  of  any  side  of 
a  plane  triangle  Is  equal  to  the  sum  of  the  squares  of  the  other  two  sides 
diminished  by  twice  the  product  of  those  two  sides  and  the  cosine  of  their 
Included  angle."  Figure  B-1  may  again  be  used  to  Illustrate  the  pertinent 
geometry  of,  and  the  arbitrary  coordinate  frame  for,  the  solution  of  the 
problem.  Figure  B-5  illustrates  the  cosine  law  technique  for  determining 
the  u-v  coordinates  of  an  aircraft  relative  to  two  sonobuoys.  The  sample 
problem  illustrating  this  technique  is  shown  in  figure  B-6. 
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Recall  that  the  cosine  law  coordinates  have  been  trasnlated  bjr  4-1,000  yd 
along  the  u  axis  relsUvs  to  the  transfornatlon  technique  coordinates. 
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T 


Law  of  Coalneai 

r|  -  rJ  +  -  2RiR5  oot  e 

fron  iA16b 


006  9 


R£  ♦  R| 


2  -  r2 

.A  **  XLa 


2R1R3 

Bia7  bo  obtained. 
Oien, 

U  ■  Rj^  008  0 

V  -  Rj  8ln  0. 


noURE  B-$  -  Qeometzy  Illustrating  the  Cosine  Law  Technique 
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V 


QtT«ni  Rj^  -  10,000  yd 
R3  -  2,000  yd 
H2  ■  8,08Ji.7  yd 


9  -  15*. 


Ihen 

u  -  Ri  008  0  -  10,000  (0.56593)  -  9,659.3 
T  -  Ri  aln  0  -  10,000  (0.25802)  -  2,588.2. 

NOTE1  In  setting  up  the  problem,  R^  and  R3  were  fixed  and  9  was  taken 

as  15  degrees,  u  and  v  were  oonputed  and  the  dlstanoe  formula 
was  employed  to  ooifnite  R2. 


JTOURE  B-6  -  Sample  Problem  Illustrating  the  Cosine  Law  Teohnlque 
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ADDENDUM 


By  a  simple  mathematloal  manipulation  it  may  be  shoun  that  Ag  and  Ay 

may  be  computed  directly  in  terms  of  R2,  and  Ro  by  the  use  of  the 
following  formula. 


Rl  +  R2 


;  Ay 


R3 
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APPENDIX  C 

MATHEMATICS  OF  TRANSFORMATION  TECHNIQUE  APPLIED  TO 
INTERSECnNG  ELLIPSES 

(The  equations  for  the  second  practicable  application  of  the  "Transforma¬ 
tion  Technique"  are  developed  herein) 
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MATHEMATICS  OP  TRANSFORMATION 
TECHNIQUE  APPLIED  TO 
INTERSECTING  ELLIPSES 


DISCUSSION 

The  problem  geometry  la  ehoim  In  figure  C>1. 


V 


FIGURE  C-1  -  Problem  Geometry 

Ihe  target  location  (u,  v)  may  be  determined  from  the  elmultaneoua 


solution  oft 

<^(u  -  b)*  ♦  (v  -  c  f  ♦  ^ (u  +  a)*  ♦  -  dj^  (C-l) 

^(u  -  b)"  +  (v  -  c)*  +  ^(u  -  a)"  +  -  dp  (C-2) 

To  solve I 

Rewrite  equation  (C-l)  as 

^(u  -  b)*  ♦  (v  -  c)*  -  dj^  -  -  ^(u  +  a)*  +  v* .  (C-3) 

Square  equation  (C-3)  to  obtain 

-2bu  ♦  I?  -  2cv  +  -  2di  ^(u  -  b)*  +  (v  -  o)* 

♦  d^  -  2au  -  «?  -  0.  (C-lt) 

Rewrite  equation  (C-2)  as 

-^(u  -  b)*  +  {y  -  cf  -  dp  "  -  ^(u  -  a)*  *  .  (C-5) 
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Square  equation  (C-5)  to  obtain 

-bftt  ♦  -  2o!ir  +  <^  -  2d2  »J  (u  -  b)*  ♦  (v  -  e)* 

♦  ^  *  200.  -  1^.  •  Oi  (C-6) 

Arrange  texlaa  In  equatlai  (C-U)  and  (6-6)  to  obtain 

(^2b  -  2a)u  -  2ev  ♦  I?  ♦  c?  ♦  d^  -V  -  2di  y^(u  -  b)®  +  (v  -  c)* 

,  .  0  (C-7) 

(-2b  +  2a)u  -  2ov  ♦  t?  *  +  d|  -  -  2d2  ^(u  -  b)*  +  (v  -  c)® 

-  0.  (C-8) 

Now,  multiply  equation  (C-7)  by  d2,  equation  (C-8)  by  d^,  pubtract 
and  arrange  terms  to  obtain 

u  [d2  (-2b  -  2a)  -  di-(-2b  ♦  2a)]  ♦  v  [-2cd2  -  (-2cdi)] 

♦  d2  (I?  ♦  ♦  df  -  rf*)  -  di  ♦  e?  +  d^  -V)  *  0.  (0-9) 

Let  R  ■  2  [b(d]^  -  d2)  -  a  (d2  -  dj^)] 

P  -  2o  [di  -  d2] 

M  -  V  ♦  cP  -  ^ 

Q  -  d2  (M  ♦  dj)  -  d^  (M  ♦  d|), 

M  obtain 

Ru  ♦  Pv  ♦  Q  ■  0.  (C-10) 

Now,  to  obtain  a  one-to-one  mapping  for  the  hypertela,  the  problem 
met  be  normalised  so  that  a  s  1.  Therefore,  the  geometry  of  flpire 
C-1  assumes  the  form  of  figure  C-2. 

From  figure  C-2  It  oan  be  seen  that 

p''  *  R^  defines  an  ellipse  -  ellipse  NO.  1 

p*  *  R2  defines  an  ellipse  -  ellipse  No.  2 

defines  a  hyperbola 
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nOURE  C-2  -  Normalized  Problem  Geometry 
These  oonlo  sections  are  Illustrated  in  figure  C-3 


V 


figure  C-3  -  Conic  Sections  Pertinent  to  Problem  Solution 
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The  transformation 
u'  ■  a^  cosh  y 
v'  -  bg  sinh  y 

where  a^  and  srQ  defined  by 
(»')*  (r'f  _  j 

(•i)* 

maps  the  rij^t  hand  branch  of  the  hyperbola,  shown  in  figure  C-3f  into 
the  Btraii^t  line  H|{  in  the  x,  y  plane  as  shown  in  figure  C-li. 


I' 

1 


X 


Image  of  point 
Image  of  point  I2 


FIGURE  C-U  -  Map  of  Right-Hand  Branch  of  the  Hyperbola  of  Figure  C-3 
into  the  Straij^t  Line  Hj^ 

NOTEi  It  is  not,  in  general,  necessary  to  know  how  the  ellipses  of  fig¬ 
ure  C-3  nap  into  the  x-y  plane  of  figure  C-U|  however,  the  points 
of  intersection  {li  and  I2)  of  ellipses  1  and  2  transform  in  the 
same  manner  as  ipoints  on  the  hyperbola.  The  inverse  mapping  of 
points  Ix  snd  I2  from  figure  C-U  will  uniquely  determine  points 
Ix  and  I2  respectively,  of  figure  C-3>  These  points  (Ix  and  I2) 
are,  of  course,  the  solutions  that  are  souf^t  after. 

Referring  to  figure  C-2,  P,  Q,  and  R  become,  respectively, 

P.2c' 

Q  -  d^  [  M  +  (di)“  1  -  dx  [m  +  (d^)*] 

R  -  2  [b'  (dx  -  dg^)  -  <4  -  dx] 
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where, 

o'  =  c/a 
b’  =  b/a 
d|  =  d^^/a 
dg  =  dg/a 

M  =  (b>)®  -  1 

and  equation  (C-IO)  becomes 

Ru'  +  Pv'  +  Q  =  0.  (C-11) 

Now,  using  the  information  implicit  in  figures  C-2  and  C-3,  replace 
u'  and  v'  in  equation  (C-ll)  by  aw  cosh  y  and  bj^  sinh  y,  respectively,  to 
obtain ; 

R  an  cosh  y  +  P  b^  sinh  y  +  Q  =  0.  (C-12) 

Then, 

R  a^  cosh  y  ±Pb|^  ^ oosh*  y  -  1  +  Q  =  0 

and 

±P  b^  .^ycosh^  y  -  1  =  -  j^R  an  cosh  y  +  Qj 
square  to  obtain 

cosh?  y  -  ■  i^(a^)®  cosh?  y  +  2QRa^  cosh  y  +  Qp 

finally, 

-  f?  cosh?  y  -  2  QRa^  cosh  y  - 

=  0.  (C-13) 

Now,  let  Z  =  cosh  y 

A  =  I^(bj^)*  -  R?  (ajJ)« 

B  =  -2QRa{J 
C  =  (bji)»  -  (f  . 
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To  obtain 

k  t  +  BZ  +  C  =  0 


where 


(C-lli) 


and 

2  °  ®H  2 
I  -Run  2  ~  ^ 

V,  o  =  - - • 

i,  p 

Finally, 

'll,  2  =  a  'll,  2 
^1,  2  =  3  '^l,  2 

It  may  be  noted  that  the  above  technique  does  not  necessitate  the 
direct  computation  of  either  a  hyperbolic  function  or  argument  thereof 
per  se. 

It  the  target  is  on  the  v-axis,  a  value  of  A  =  0  will  result  and  the 
expression  for  Z  becomes  infinite.  Tb  avoid  this  complication,  the  fol¬ 
lowing  alternate  must  be  used. 


If  A 

then 


Hi,  2 


0 

=  0 


where  p 

G  =  ii  (c  -  d£) 

I 

di  *  ad^  or  +  a 

M  =  t?+C*-8(®-di 
H  -  -li  c  M 
Z  =  M®  -  lidjaf* 

Ihe  subscript  1  here  implies  use  of  the  (+)  sign. 
The  subscript  2  here  implies  use  of  the  (-)  sign. 
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APPENDIX  D 


PROGRAMMING  ANALYSIS 
OF 

APPLICATIONS  I  AND  II 
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PROGRAMMING  ANALYSIS  OF 
APPLICATIONS  I  AND  II 

Submitted  by  E.  A.  Book,  Jr 

f 

Since  there  may  be  more  than  one  acceptable  method  of  solution  avail¬ 
able  (from  a  strictly  m,athematical  point  of  view)  for  a  given  problem, 
it  was  necessary  to  utilize  a  methodology  for  choosing  that  solution 
which  would  be  the  most  efficient  when  machine  computation  is  to  be  con¬ 
sidered.  The  criteria  used  were  the  following: 

1.  Precision  of  the  solution  (based  on  a  21  bit  data  word) 

2.  Program  running  time 

3 .  Space  required 

il.  Error  analysis 

The  final  program  of  each  method  was  run  on  a  PB250  computer  to  verify 
the  results  of  the  methodology. 

The  final  analysis  showed  that  the  transformation  technique  was  the 
most  efficient  method  to  use  with  a  computer.  Less  program  analysis 
was  required  with  the  transformation  techjiiqae  than  with  either  of  the 
other  two  methods  tried.  Tables  I>-I  and  D-II  show  the  results  of  the 
evaluation,  and  clearly  indicate  that  the  transformation  technique  is 
superior  to  the  other  techniques. 

TABLE  D-I 


Type  of 

Accuracy 

Running  Time 

Space 

Solution 

(Decimal  Digits) 

(usee) 

(21  bit  words) 

Cosine  Law 

k 

3,21;0 

57 

Transforma  tion 

5 

3,U56 

56 

TABLE 

D-II 

Type  of 

Accuracy 

Running  Time 

Space 

Solution 

(Decimal  Digits) 

(usee) 

(21  bit  words) 

Ellipse-Ellipse 

3 

11,256 

192 

Transformation 

3 

9,528 

130 
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DETAILED  RESULTS  OF  THE  ANALYSIS  OP  mE  COSINE  LAW  AND  TRANSFORMA  HON 
TECHNIQUES 

The  nvuiiber  of  Instructions  and  the  running  time  are  essentially  the 
same  with  both  techniques.  However,  more  fundamental  analysis  was 
required  with  the  cosine  law  technique,  because  an  Increased  amount  of 
scaling  was  necessary  to  produce  an  optimum  solution. 

The  analysis  used  to  determine  the  maximum  number  of  significant 
decimal  digits,  which  would  be  obtained  from  each  solution,  revealed 
that  the  transformation  technique  result  contained  at  least  one  more 
significant  dlgltj  therefore.  It  would  produce  more  accurate  results 
(with  the  same  Input)  than  the  cosine  law  technique. 

DETAILED  RESULTS  OF  THE  ANALYSIS  OF  THE  ELLIPSE-ELLIPSE  AND  TRANSFOR¬ 
MATION  TECHNIQUES 

The  analysis  showed  that  even  when  the  Input  data  was  normalized, 
the  competitive  technique  required  a  "floating  point"  routine  to  produce 
a  solution.  Although  ^e  use  of  the  floating  point  routine  reduced  the 
analysis  required  for  scaling.  It  slowed  down  the  running  time  needed  to 
produce  a  solution.  This  Increase  In  time  Is  not  reflected  In  table 
D-II;  this  was  not  the  case  with  the  transformation  technique. 

The  running  time  for  the  ellipse-ellipse  technique  was  much  longer 
than  for  the  transformation  technique.  This  was  due  partly  to  the 
Increased  number  of  Instructions  needed  and  partly  to  an  Increased 
number  of  divide,  multiply  and  square  root  Instructions  needed. 

With  the  criteria  used,  the  resvilts  of  the  analysis  show  that  the 
transformation  technique  was  the  most  efficient  program  to  use  when 
solving  for  the  points  of  Intersection  of*  conic  sections  more  complex 
than  circles. 
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